Mesoscopic supercurrent transistor controlled by nonequilibrium cooling 
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The distinctive quasiparticle distribution existing under nonequilibrium in a superconductor- 
insulator-normal metal-insulator-superconductor (SINIS) mesoscopic line is proposed as a novel 
tool to control the supercurrent intensity in a long Josephson weak link. We present a description 
of this system in the framework of the diffusive-limit quasiclassical Green-function theory and take 
into account the effects of inelastic scattering with arbitrary strength. Supercurrent enhancement 
and suppression, including a marked transition to a 7r-junction are striking features leading to a 
fully tunable structure. The role of the degree of nonequilibrium, temperature, and materials choice 
as well as features like noise, switching time, and current and power gain are also addressed. 

PACS numbers: 74.50.+r, 73.23.-b, 74.40.+k 
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I. INTRODUCTION 



Nonequilibrium effects in mesoscopic superconduct- 
ing circuits have been receiving a rekindled attention 
during the last few years 0. The art of control- 
ling Josephson coupling in superconductor-normal metal- 
superconductor (SNS) weak links is at present in the 
spotlight: a recent breakthrough in mesoscopic super- 
conductivity is indeed represented by the SNS transistor, 
where supercurrent suppression as well as its sign rever- 
sal (7r-transition) were demonstrated 0, 0, 0, 01 Hi . This 
was achieved by driving the quasiparticle distribution in 
the weak link far from equilibrium [3j ISl l3 through ex- 
ternal voltage terminals, viz. normal reservoirs. Such a 
behavior relies on the two-step shape of the quasiparti- 
cle nonequilibrium distribution, typical of diffusive meso- 
scopic wires and experimentally observed by Pothier and 
coworkers [Tof . 

The purpose of this paper is to demonstrate that it is 
possible to tailor the quasiparticle distribution through 
superconductivity-induced nonequilibrium in order to im- 
plement a unique class of superconducting transistors 
hll. This can be achieved when mesoscopic control lines 
are connected to superconducting reservoirs through tun- 
nel barriers (I), realizing a SINIS channel. The pecu- 
liar quasiparticle distribution in the N region, originating 
from biasing the S terminals, allows one to access several 
regimes, from supercurrent enhancement with respect to 
equilibrium to a large amplitude of the 7r-transition pass- 
ing through a steep supercurrent suppression. These fea- 
tures are accompanied by a large current gain (up to some 
10 5 in the region of larger input impedance) and reduced 
dissipation. The ultimate operating frequencies avail- 
able open the way to the exploitation of this scheme for 
the implementation of ultrafast cryogenic current and/or 
power amplifiers. 

The paper is organized as follows. In Sec. II we in- 
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troducc diffusive-limit quasiclassical Green-function the- 
ory that is employed to describe the supercurrent in the 
proposed structure. In particular, we show how to ex- 
ploit superconductivity-induced nonequilibrium in order 
to control the Josephson current. In Sec. Ill we address 
the role of inelastic scattering in the structure by solv- 
ing the kinetic equation for the SINIS control line, and 
we describe how the strength of such scattering affects 
the observable supercurrent. The critical current depen- 
dence on the length of the weak link is analyzed in Sec. 
IV where the relevance of long junctions is pointed out. 
In Sec. V we consider the effect of the additional control 
terminals on the observable supercurrent while Sec. VI 
is devoted to the analysis of the supercurrent harmonics. 
In Sec. VII we address the role of materials combina- 
tion in determining the performance of the structure. In 
particular, issues like power dissipation and noise power 
are discussed. The achievable current and power gain 
as well as the maximum attainable transistor operating 
frequency are discussed in Sec. VIII. Finally, Sec. IX 
summarizes the main results. 



II. TRANSISTOR OUTPUT CHARACTERISTIC 

The investigated mesoscopic structure (see Fig.^) con- 
sists of a long diffusive weak link of length Lj much larger 
than the superconducting coherence length (£o) oriented 
along the x direction. This defines the SNS junction of 
cross-section Aj. The superconducting terminals belong- 
ing to the SNS junction, labelled Sj (3 and 4), are kept at 
zero potential in the analysis of the supercurrent (transis- 
tor operation may take place in the dissipative regime at 
a finite voltage) . The SINIS control line is oriented along 
the y direction and consists of a normal wire, of length Lc 
and cross-section Ac , connected through identical tunnel 
junctions of resistance TZt to two superconducting reser- 
voirs Sc (1 and 2), biased at opposite voltages ±Vc/2. 
The superconducting gaps of Sj and Sc (Aj and Ac) 
are in general different in order to optimize the structure 
characteristics. 

The supercurrent Ij flowing across the SNS junction 
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FIG. 1: Scheme of the Josephson transistor. The supercurrent 
I j (along the white dashed line) is tuned by applying a bias 
Vc across the SINIS symmetric line connected to the center 
of the weak link. The superconducting gaps Aj and Ac are, 
in general, different, and all normal wires are assumed quasi- 
one- dimensional, i.e., their width is much smaller than their 
length. 

is given by d, || 

(7 A f°° 

h{Vc) = -r 1 dE [f(-E; V c ) - f(E; V c )]Im\j E ], 

(1) 

and depends on the quasiparticle distribution function 
f(E) in N. In Eq. (JIJ, a is the normal-state conductiv- 
ity which determines the normal-state resistance of the 
junction according to = Lj/aAj. As shown below, 
the distribution function / is determined by the SINIS 
control line and reduces to the equilibrium Fermi distri- 
bution when Vc = 0. 

The energy-dependent spectral supercurrent, 0, H^ | 
Imb'rej, can be calculated by solving the Usadel equa- 
tions [14) . Following the parametrization of the Green 
functions given in Ref. 0, these equations in the N re- 
gion can be written as 

j E = -sinh 2 (6)3 xX , d x j E = 0, (2) 
hDd 2 x 9 + 2iEsmh9 - ^-(d xX ) 2 sinh(2<9) = 0, (3) 

where D — \ypi m is the diffusion coefficient, Vp is the 
Fermi velocity, £ m is the mean free path and E is the 
energy relative to the chemical potential in Sj. 9{x,E) 
and X (x,E) are in general complex functions. For per- 
fectly transmissive contacts, the boundary conditions at 
the SjN interfaces reduce to 9 = arctanh(Aj/.E7) and 
X = ±0/2 in the reservoirs Sj, where (j) is the phase 
difference between the superconductors. The presence 
of the control wires on the spectral supercurrent can 
be taken into account explicitly including two additional 
probes and imposing boundary conditions at the ScN in- 
terfaces as in Ref. [l3|] ■ In the present case such boundary 
conditions take the form d x 6 = to describe large tunnel 
barriers of resistance TZt- The large TZt allows us to ne- 
glect the tiny oscillating supercurrents flowing between 




FIG. 2: (a) Supercurrent vs control voltage Vc at different 
temperatures (T) for <f> = tt/2 and T c c = 0.2 T/ (see text). 
Bias regions labelled (I), (II) and (III) indicate supercur- 
rent enhancement due to quasiparticle cooling, high- voltage n- 
state and low-voltage 7r-state in the high-temperature regime, 
respectively. These are qualitatively explained in (b), (c) and 
(d) where hatched areas represent the contribution to super- 
current arising in such bias ranges (see text). We assume 
throughout that supercurrent at <j> = n/2 yields values close 
enough to the critical current, an assumption that breaks 
down near the ir -junction transition [(J (see Sec. IVH . 



the superconductors Sc- As discussed in Sec. V, the su- 
percurrent Ij depends trivially on the control wire length 
Lc, namely it monotonically decreases, rapidly reaching 
an asymptotic value of, roughly, one half the zero-length 
value. This in turn allows us to make the following sim- 
plifications: as far as the spectral current is concerned, 
we shall assume Lc/Lj — * 0, aware of the fact that the 
effect of a finite Lc/Lj gives rise to a mere reduction of 
Ij. A similar discussion can be made on the dependence 
of the equilibrium supercurrent on Ac- As shown in Ref. 
[l3| , Ij is a monotonically decreasing function of Ac /Aj 
and, as far as spectral supercurrent is concerned, we as- 
sume Ac/Aj — > throughout the paper (in fact, it is 
enough to assume either Lc *C Lj or Ac -C Aj). 

We must now determine the actual quasiparticle distri- 
bution / in the N region of the SINIS structure. This is 
controlled by voltage (Vc) and temperature, and by the 
amount of inelastic scattering in the control line. In the 
case of a short control wire with no inelastic interactions, 
the quasiparticle distribution, according to Ref. [15) . is 
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riven by 



III. ROLE OF INELASTIC SCATTERING 



f(E,V c ) 



M + A/a 



(4) 



where A/i, 2 = N Sc ( E ± eVc/2) and JF 1;2 = JF°(£7 ± 
eVc/2). The former are the BCS density of states in the 
reservoirs Sc (labelled 1 and 2 in Fig. QJ. F°(E) is the 
Fermi function at lattice temperature T 01- Note that 
expression |(2J can also be found from the quasiclassical 
theory in the tunnelling limit [lfj. In this case Eqs. 
and (QJ yield the dimcnsionlcss transistor output charac- 
teristics shown in Fig. |2fa), where we plot the supercur- 
rent Ij versus control bias Vc at different temperatures. 
We assumed 4> = tt/2, T c c /T c j = 0.2, where T C C(J) are 
the critical temperatures of the superconductors Sc(j), 
and Lj such that Aj/E^ h = 300. We choose the limit 
of a long junction (i.e., Aj ^> where E^ = HD/Lj 
is the Thouless energy of the SNS junction), since in this 
limit the supercurrent spectrum varies strongly with en- 
ergy. 

For all temperatures T < TjP, with increasing Vc, 
curves display a large supercurrent enhancement with re- 
spect to equilibrium at bias V c = 2A c {T)/e = V C (T) 
(region I in the figure). Further increase of the bias 
leads to a 7r-transition (region II) and finally to a de- 
cay for larger voltages [18|- This behavior is illustrated 
by Figs. 2(b,c,d) where the spectral supercurrent (solid 
line) is plotted together with f(—E) — f(E) (dash-dotted 
line) for values of Vc and T corresponding to regions I, 
II and III, respectively. Hatched areas represent the in- 
tegral of their product, i.e., the supercurrent Ij of Eq. 
(|T|). In particular, region I corresponds to the cooling 
regime where hot quasiparticles are extracted from the 
normal metal [TH |20| . The origin of the 7r-transition in 
region II is illustrated by Fig. Hfc), where the negative 
contribution to the integral is shown. We remark that 
the intensity of the supercurrent inversion is very sig- 
nificant. It reaches about 60% of the maximum value 
of Ij at Vc — V C (T) in the whole temperature range, 
nearly doubling the 7r-state value of the supercurrent as 
compared to the case of an all-normal control channel 
[llg. In the high-temperature regime (T/Xf > 0.6), 
when the equilibrium critical current is vanishing, the 
supercurrent first undergoes a low-bias 7r-transition (re- 
gion III in the figure), then enters regions I and II. This 
recovery of the supercurrent from vanishingly small val- 
ues at equilibrium is again a consequence of the peculiar 
shape of / (see Fig. EJd)). Notably, the supercurrent 
enhancement around V C (T) remains pronounced even 
at the highest temperatures, so that Ij attains values 
largely exceeding 50% of the junction maximum super- 
current. This demonstrates the strong tunability of the 
supercurrent through nonequilibrium effects induced by 
the superconducting control lines. We remark that this is 
a unique feature stemming from the superconductivity- 
induced nonequilibrium population in the weak link. 



Since we do not want to confine our analysis only to the 
limit of control wires much shorter than the energy relax- 
ation length, we need to evaluate the impact of inelastic 
scatterin g pr esent in the N region on the supercurrent 
behavior |2l|. As a matter of fact, the length Lc of the 
SINIS control line can be additionally varied to control 
the supercurrent by changing the effective strength of in- 
elastic scattering in the N region. In general, the steady 
state distribution function at a given energy E is deter- 
mined by an equation which accounts for the balancing of 
quasiparticle flow to and from the reservoirs with a colli- 
sion term I, describing energy relaxation processes. For 
TZt ^S> Rc = Lc/crAc, the distribution function f(E) in 
the N region is essentially y-independent and we have 
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e 2 TZT^lc v F 



[M(^i - f{E))+M 2 {T 2 - f(E))] + 
+k J dudeuj a l{uj,£,E) = 0. 



(5) 



Here vf is the normal-metal density of states at the Fermi 
energy, £lc is the volume of the N region and X is the 
net collision rate at energy E. At low temperatures, the 
most relevant scattering mechanism is electron-electron 
scattering j2^ and we can ne glec t the effect of electron- 
phonon scattering. Then, [Toll23| 



l{uj,e,E)=l m {Lu,e,E)-l out {uj,e,E), 



(6) 



and 



T m (iu, e, E) = (l- f(e))(l - f(E))f(e - u)f(E + W ),(7) 
I mit {u, s, E) = (1 - f(e - w))(l - f(E + w))/( e )/(#).(8) 

Electron-electron interaction is either due to direct 
Coulomb scattering |2J, |25j or mediated by magnetic 
impurities |22j . Below, we concentrate on the former, 
but the latter would yield a similar qualitative behavior. 
From the calculation of the screened Coulomb interaction 
in the diffusive channel, it follows [24| that a = —3/2 for a 
quasi-one dimensional wire and k = {tx\J2DY?I 2 vfAc)~ 1 
[25l |2(| . We note that Ac is the most relevant energy 
scale to describe the distribution function for different 
voltages Vc- It is thus useful to replace u — > uj/Ac and 
e — * ej Ac, in order to obtain a dimensionless equation. 
Multiplying Eq. 10 by e 2 Tl T ?LcVF we obtain 



M X (T X - f(E)) - N 2 {f{E) - T 2 ) = 
= fC co ii J du;deuj~ 3/2 l{Lj,£,E), 



where 




(9) 



(10) 



Rk = h/2e 2 and E^ h — hD/L c . In the absence of 
electron-electron interaction (Kcoll = 0), Eq. Q is re- 
covered. 
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FIG. 3: Supercurrent vs Vc for various /C co ;; with T = 0.1 T,r 
and = 0.2 T C J . Insets show the distribution function at 
eVc = 1.5A C (left) and eVc = 2.5A C (right) calculated for 
the same K. co u values. 



FIG. 4: Effect of different ratios Aj/E^ h on the observable 
supercurrent at T = 0.1 T c c and T c c = 0.2 T/. Note that the 
effect of the nonequilibrium distribution on the observable 
supercurrent is more pronounced in longer junctions. 



The influence of inelastic scattering on Ij is shown 
in Fig. |21 which displays the critical current of a long 
junction at T = 0.1 T c c for several values of K, C oii- Here Ij 
is obtained by numerically solving Eq. © . The effect of 
electron-electron interaction is to strongly suppress the ir- 
state and to widen the peak around V£ . The 7r-transition 
vanishes for K. co n ~ 100, but the Ij enhancement due 
to quasiparticlc cooling still persists in the limit of even 
larger inelastic scattering [23 ■ 

The disappearance of the 7r-state can be understood by 
looking at the right inset of Fig. [21 which clearly shows 
how a large-bias / (calculated at eVc = 2.5A C ) gradu- 
ally relaxes from nonequilibrium towards a Fermi func- 
tion upon increasing K co u- The left inset shows how a 
low-bias / (evaluated at eVc — 1.5A C ) sharpens, thus 
enhancing Ij, by increasing K, co ii- This effect follows 
from the fact that inelastic interactions redistribute the 
occupation of quasiparticlc levels in the N region, thus 
increasing the occupation at higher energy. As a con- 
sequence, higher-energy excitations are more efficiently 
removed by tunneling, even for biases well below and not 
only around V£ (as in the case of K, co u =0). At the 
same time, supercurrent recovery at high temperature 
is gradually weakened upon enhancing IC co u. Notably, 
these calculations show that a rather large amount of in- 
elastic scattering is necessary to weaken and completely 
suppress the 7r-state. For example, using AI/AI2O3/CU 
as materials composing the SINIS line, K co ll = 1 corre- 
sponds to use a fairly long control line with Lq — 4.6 /im 



IV. DEPENDENCE ON JUNCTION LENGTH 

The choice of the superconductors Sj and Sc affects 
the behavior of the device through the energy scales given 
by the energy gaps Aj/c in these materials and their ra- 



tio with the geometric scales of the system. The role 
of junction length [Lj) on the transistor characteristics 
is displayed in Fig. ^ where the supercurrent is plot- 
ted for different ratios Aj/E^ h at T = 0.1 T c c and for 
T,r = 0.2 T/. Upon shortening the junction the Joseph- 
son current response to the control voltage Vc becomes 
wider (note that also the ratio between the characteristic 
energies Ac and Ei, h changes — our model holds only for 
A j > Ac). At the same time, features like supercurrent 
enhancement and 7r-transition appear less pronounced. 
This suggests that large Aj/E^ h ratios are favorable for 
an efficient transistor effect, so that the Ij(Vc) charac- 
teristic is forced into a narrower bias window. 



V. EFFECT OF EXTRA TERMINALS 

As mentioned earlier, the presence of the lateral arms 
influences the superconducting correlations induced in 
the weak link by reducing 13] the spectral supercurrent. 
Such an effect is shown in Fig. where the equilibrium 
supercurrent (i.e., at Vc = 0) is plotted as a function 
of the ratio Lq / Lj. The magnitude of the supercurrent 
is decreased upon increasing Lc, reaching its asymptotic 
value at Lc — Lj, almost independently of temperature. 
Note that, for each temperature, such a value is roughly 
one half that corresponding to the limit Lc — > 0. This 
conclusion contrasts with the case of a normal control 
line without tunnel barriers 0, ^| and emphasizes the 
relevance of the condition Lc *C Lj towards the maxi- 
mization of the supercurrent in the system. 



VI. PHASE DEPENDENCE 

The generic sinusoidal supercurrent-phase dependence 
is strictly valid only in the case when the Josephson junc- 
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FIG. 5: Equilibrium supercurrent vs Lc/Lj at different tem- 
peratures. Note that 7j reaches its asymptotic value at 
Lc — L.j, almost independently of the temperature. 



FIG. 6: (a) Voltage dependence of the amplitude of the first 
two harmonics of the observable supercurrent calculated at 
T = O.lTf and Tf = 0.2 T C J . (b) The same as in (a) calcu- 
lated in the high-temperature regime (T = 0.8 T c c ). 



tion is a tunnel junction. Transmission channels with 
high transmission probability, present in diffusive wires, 
tend to deform this dependence by bringing in higher 
harmonics. Thus, the supercurrent-phase relation is in 
general of the form [l3| 

oo 

W) = X^,» sin N>), (ii) 

71=1 

where are the coefficients of the Fourier sine series of 
Is{4>)- Loosely speaking, the higher harmonics arise from 
multiple passages of Cooper pairs through the junction 
in a single coherent process. In most cases for a diffu- 
sive junction, the n = 1 harmonic is dominant, and the 
role of the others is to slightly deform the current-phase 
relation, without bringing about additional nodes. How- 
ever, across the 7r-junction transition, the first harmonic 
can be completely suppressed. In this case, the current- 
phase relation is determined by the second harmonic, and 
around this point, the supercurrent has a half periodicity 
with respect to the phase (flux in a loop). Such a sin(2</>) 
phase dependence was measured in Ref. jg in the case 
of normal-metallic control probes. 

Figure [S] shows the voltage dependence of the two low- 
est harmonics at two temperatures. One can find that 
for most voltages, the first harmonic is dominant, but in 
the region(s) of the transition between the conventional 
and the 7r-states, the second becomes dominant. More- 
over, one finds that, analogously to the case with normal 
control probes, 0, ^1 the S1 g n of the second harmonic 
in these regions is positive. Hence, in those values of the 
voltage Vc, the free energy of thejunction has local min- 
ima both at <f> = and <j> = tt |6(|. Another implication 
of the finite value of Ij t 2 around those voltages is that 
the critical current of the junction does not completely 
vanish, but it is merely no longer obtained for <fi sa tt/2. 



VII. PERFORMANCE I: ROLE OF MATERIALS 
CHOICE 

It is interesting to analyze the role of different su- 
perconductor combinations, i.e., different ratios X^/T/. 
This effect is shown in Fig. where the supercurrent vs 
bias voltage (normalized to Aj/e) characteristic is dis- 
played for different combinations of superconductors at 
T = 0.02 T/. The main effect is a shift in the Ij response 
along the Vc axis towards lower bias voltages upon re- 
ducing the critical temperature T c c . Furthermore, the 
shape of the characteristics appears to be virtually inde- 
pendent of the TjF value, although a small decrease of the 
7r-transition amplitude gradually develops by decreasing 
the ratio T^/T/. Such a ratio can be chosen to move 
the range of relevant Vc towards lower values and thus 
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FIG. 7: Supercurrent vs Vc calculated at T = 0.02 T C J for 
different Tf/T/ ratios. The main effect of lowering is to 
shift the I j response along the Vc axis and to slightly reduce 
the amplitude of the 7r-transition. 
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FIG. 8: (a) Power dissipated in the SINIS line vs Vc calcu- 
lated for various ratios Tf /T C J and T = 0.01T/. (b) Noise 
power S vs Vc calculated for the same parameters as in (a). 
In all these calculations we set T C J = 9.26 K corresponding to 
Nb and 11t = \ kQ. 



FIG. 9: Differential current gain G/ vs Vc for T c c /T/ = 0.2 
and T = 0.01 T C J : in (a) G/ is shown for Vc < VS while in 
(b) it is shown in the high-bias region. In these calculations 
we set Kcoll = and T C J = 9.26 K (Nb). 



decrease the power dissipation P = IcVc, where Ic is 
the control current across the SINIS channel. The func- 
tion P(V C ) is plotted in Fig.^a) for some ratios T c c /T/ 
at T — 0.01 T/, assuming TZ T = 10 3 Q and T C J = 9.26 
K (corresponding to Nb). The impact of Ac in control- 
ling power dissipation is easily recognized. These effects 
clearly point out that Ac <§C Aj is the condition to be 
fulfilled in order to minimize P. In practice, the power 
dissipation for eVc > 2A constitutes an experimental 
problem as this energy needs to be carried out from the 
reservoirs. 

In a similar way the noise properties of the system are 
sensitive to the different X^/T/ ratios. Assuming that 
the noise through one junction is essentially uncorrelated 
from the noise through the other, it follows that the input 
noise power S in the control line can be expressed as 

i r°° 

S(V C ) = — / dEM 1 [f{E){l-T 1 )+T 1 {l -f(E))}. 

(12) 

The function S(Vc) is shown in Fig. Elb) for the same 
parameters as in Fig. |Sfa). For example, for T^/T/ 
0.1 (corresponding roughly to the combination Al/Nb), 
P obtains values of the order of a few fW and S of some 
10~ 30 A 2 Hz _1 in the cooling regime, while these values 
are enhanced respectively to few tens of pW and 10~ 26 
A 2 Hz _1 for biases around the 7r-transition. 



Af/e and Ej, h /e, respectively, so that the voltage gain 
is G v = V out /V in ra E J T JA C . As A c -> E J Thl the de- 
vice characteristics become smoother (see Sec. IIV() and 
the control becomes less efficient, so that it is hardly 
possible to achieve voltage gain. On the other hand, dif- 
ferential current gain G/ = dlj/dlc can be very large. 
For Vc > Vq, a simple estimate gives 

Gi ~ (E^/AcKTIt/Rn), (13) 

meaning that with realistic ratios TZt/Rn (~ 10 3 ), Gj 
can exceed 10 2 . G](Vc) can be calculated noting that 
Gj = (^)C^)- 1 . In Fig.Hwe plot G, for Tf/T/ = 
0.2. This calculation reveals that G r ca n reach huge val- 
ues, with some 10 s for Vc < V£ and several 10 2 
in the opposite regime. Remarkably, gain is almost un- 
changed also in the presence of weak inelastic scattering 
(i.e., JCcoii = 1). The same holds for P and S. 

Operation of the Josephson junction in the dissipa- 
tive regime makes it possible to exploit the system for 
power gain. Let us assume that the voltage across the 
Josephson junction (Vj) is fixed and that the current- 
voltage characteristic is nonhysteretic. According to the 
RSJ (resistively shunted junction) model [3(j for over- 
damped junctions, the voltage drop across the junction 
at T — reads 



VIII. PERFORMANCE II: GAIN AND 
SWITCHING TIME 

In light of the possible use of this operational prin- 
ciple for device implementation, let us comment on the 
available gain and switching times. Input (Vi„ = Vc) 
and output (V ou t — IjRn) voltages are of the order of 



Vj = R Ny Jl*-lj c , (14) 

where Ij c is the critical current and / is the total current 
through the Josephson junction. The differential power 
gain is defined as Gp = dP ou t/dPi n where P ou t — Vjl, 
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Pin — Vole- Using Eq. ifH|) one obtains 

dl Vjlj c dI Jc dl c 



G P = Vj 



d(V c I 



c 



\Rn ) 



a ~ die d(V c I c ) 

+ l2 Jc 



'Jc 



Vj 



2 „ Vc + IcRswis(Vc 



(15) 

-Gj, (16) 



where i?siNis = dVc / die is the differential resistance of 
the SINIS line at the operating current Iq- Choosing 
Vj = IjcRn we obtain 



Vj 



Gj > 



Vj 



V2(V c +I c Rsmis) ~ 2V2Vc 



(17) 



The thermal smearing of Eq. I|14|) 30] will somewhat de- 
crease Gp, but not by orders of magnitude. A further 
simplification of Eq. (|17f) can be additionally given, re- 
calling that Vj ~ E'Thl e an d Vc ~ 2Ac/e, so that 



4V2A ( 



-G/ ~ GyG/. 



(18) 



A straightforward estimate for the differential power gain 
from Eq. JTSJl gives G P - 10 2 10 3 for V c < V£ and 
~ 10 for Vc > Vq. This means that, in the suitable 
operating mode, the present structure is able to provide 
large power gain. 

As a final issue we want to briefly address the 
time scales intrinsic to this mesoscopic device when 
it is operated in the dynamic regime. The high- 
est operating frequency v of the transistor is lim- 
ited by the smallest energy in the system: v < 
mini{A c ,Aj, E^E^^UtC)- 1 ,he/ I c } where C is 
the tunnel junction capacitance. For an optimized device, 
working frequencies of the order of 10 11 Hz can be exper- 
imentally achieved in the high- voltage regime Vc > Vq. 
For Vc < Vq, conversely, the response is slower (of the 
order of some 10 8 Hz), owing to the large value of subgap 
resistance that becomes the speed-limiting factor due to 
the long discharging time through the junctions. How- 
ever, this analysis relies on the fact that the relaxation 



time scales in the superconducting reservoirs in the con- 
trol part of the structure are much shorter than the above 
scales. These scales depend on the geometry of the sys- 
tem (the real size of the reservoirs), but in practice, they 
may seriously limit the speed of the device unless special 
care is taken. We stress that the detailed prediction of 
the transistor dynamic response would require a separate 
analysis that is beyond the scope of the present paper. 



IX. CONCLUSIONS 

In this paper we have shown that superconductivity- 
induced nonequilibrium can be used to finely control the 
supercurrent flowing through a SNS Josephson junction. 
The peculiar quasiparticle energy distribution, realized in 
the normal island of a SINIS control line when a bias volt- 
age is applied between the superconducting electrodes, 
allows us to drive the SNS junction into different tran- 
sistor regimes. With increasing control voltage, we have 
found that the supercurrent first shows a steep enhance- 
ment, then undergoes a 7r-transition and finally decays 
to zero for larger voltages. Furthermore, we have ana- 
lyzed in detail the effect of inelastic scattering within the 
normal island, the dependence of the supercurrent on the 
Josephson junction length, the effect of the presence of 
lateral arms on the spectral supercurrent, and the be- 
havior of the supercurrent harmonics. Finally we have 
addressed the optimization of the device, both in terms 
of material choice and characteristic figures, such as cur- 
rent gain, power gain and operating frequencies. In view 
of these last points, we wish to point out that such a 
device could be successfully exploited as a first amplifi- 
cation stage in SQUID-based cryogenic electronics. 
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